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FOREWORD 


1. REVIEW AND PROBLEM INDEX 


Project Whirlwind makes use of the facilities of the Digital Computer 
Laboratory. Thr principal objective of the Project Is the application of an 
electronic digital computer of large capacity and very high speed (Whirlwind I) 
to problems in mathematics, science, engineering, simulation and control. 

Whirlwind I is of the high-speed electronic digital type, in which 
quantities are represented as discrete numbers and complex problems are solved by 
the tepeated use of fundamental ariiluaetic and logical (i.e., control or selection) 
operations. Computations are executed by fractional-microsecond pulses in elec¬ 
tronic circuits, of which the principal ones are (1) the flip-flop, a circuit con¬ 
taining two vacuum tubes s«* connected that one tube or the other 19 conducting, 
but not both; (2) the gate or coincidence circuit; (3) the magnetic-core memory, 
in which binary digits are stored as one of two directions of magnetic flux within 
ferro-magnetlc cores. 

Whirlwind I uses numbers of 16 binary digits (equivalent to about 5 
decimal digits). This length was selected to limit the machine to a practical 
size, but it permits the computation of many simulation problems. Calculations 
requiring greater number length are handled by the use of multiple-length numbers. 
Rapid-access magnetic-core memory has a capacity of 32,768 binary digits. Present 
speed of the computer is 40,000 single-address operations per second, equivalent 
to about 20,000 multiplications per second. 


This report covers the specific period of December 24, 1956 to March 17, 
1957. During this time, 71 problems made use of 289.51 hours of the 342.66 hours 
of Whirlwind computer time allocated to the Scientific and Engineering Computa¬ 
tions (SfcEC) Group. Of the 342.66 allocated hours of computer time, 2.57% was 
down time because of computer malfunctions. 44.35 hours of the allocated time 
were used for terminal equipment testing and calibration, demonstrations and 
various lnter-run operations not logged to specific problems. 

The 71 problems run during this quarter cover some 15 fields of appli¬ 
cation. The results of 21 of the problems have been or will be Included in aca¬ 
demic theses. In these 21 problems, there are represented 16 Doctoral theses, 3 
Master's, 2 Naval Engineer's and 1 Bachelor's. Twenty-five of the problems have 
originated from research projects sponsored at MIT by the Office of Naval Research 

Two tables are provided as an index to the problems for which progress 
reports have been submitted. In the first table, the problems are arranged ac¬ 
cording to the field of application and the source and amount of time used on 
WWI are given. In Table II, the problems are listed according to the principal 
mathematical problem Involved in each. In each table, the letter after the 
problem number indicates whether the problem is for academic credit and whether 
it is sponsored. The code is explained in section 2.1, Introduction 
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WHIRLWIND CODING AND APPLICATIONS 


2.1 Introduction 


Progress reports as submitted by the vartous programmers are presented 
in numerical order In Section 2.2. Letters have been added to the problem ‘' u “- 
bers to indicate whether the problem is for academic credit and whether it is 
SDonsozod. The letters have the following significance. 


A implies 
B implies 
C implies 
D implies 


the problem 
the problem 
the problem 
the problem 


is NOT for academic credit, is UNsponsored. 
IS for academic credit, is UNsponsored. 
is NOT for academic credit, IS sponsored. 

IS for academic credit, IS sponsored. 


126 C. DATA REDUCTION 

Problem 126 is a very large data-reduction program for use in the Servo¬ 
mechanisms Laboratory. The overall problem is composed of many component sections 
which have been developed separately and are now being combined into complete pro 
totype programs. Descriptions of the various component sections have appeared in 
past quarterly reports. After the development and testing of the prototype Whirl¬ 
wind programs is completed, the programs will be re-coded for other, commercially 
available, large scale computers, (probably the ERA 1103, IBM 701 and IBM 704 com¬ 
puters), for use by Interested agencies for actual data reduction at other loca¬ 
tions. The programs are currently being developed by Douglas T. Ross, David F. 
McAvinn, Benson H. Scheff and Dorothy A. Hamilton. Servomechanisms Laboratory 
staff members with the assistance of John F Walsh. This work Is sponsored by the 
Air Force Weapons Guidance Laboratory through DSR Project 7668. 


N implies the problem is sponsored by the Office of Naval Research. 

L Implies the problem is sponsored by Lincoln Laboratory. 

The absence of a letter Indicates that the problem originated within the S&EC 
Group. 

2.2 Progress Reports 

100 COMPREHENSIVE SYSTEM OP SERVICE ROUTINES 


The nature of the problem requires extreme automatlclty and efficiency 
in the actual running of the program, but also requires the presence of human opera 
tors in the computation loop for the purpose of decision making and program modi¬ 
fication. For this reason extensive use is made of output oscilloscopes so that 
the computer can communicate with the human, and manual intervention registers so 
that the human can communicate with the computer in terms of broad ideas, while 
the computer is running, and have the computer program translate these ideas into 
the detailed steps necessary for program modification to conform to the human 
operator's decision. The prograa> which does this translation and modification is 
called the Manual Intervention Program (MIV). The most recent version of the pro¬ 
totype data-reductlon program is called the Basic Evaluation Program. 


Staff time used on problem 100 during the past quarter has been mainly 
consumed by consultation with users of the system and with preparation of anno¬ 
tated copies of the programs involved. 

Frank Helwig 

Digital Computer Laboratory 

106 C. MIT SEISMIC PROJECT 

During the past quarter, a program for testing the hypothesis of sta- 
tionarity of a data sequence has been completed and used on several sets of data. 
Examples of random elastic media with variable scattering power have also been 
computed Finally, a spherical Bessel transform program, in WW1 code, has nearly 
been completed. 

Plans for the future include the use of the hypothesis testing program 
on a greater variety of Geophysical data and the use of the transform program to 
test the possibility that one dimensional data come from Isotropically random 
throe dimensional data. 

Programmers working on this problem are S. Simpson, S. Treltel and D. 

Grlne. 

S. Simpson 

Geology and Geophysics 


During the past quarter, work has continued in the development of SLURP 
(Servo Lab Utility Routine Program). A new generalized Group Control program Has 
been designed to take advantage of the Increased size of the magnetic core memory 
of the Whirlwind Computer. The function of the Group Control Program is to allow 
programs that are much larger than the high speed memory capacity to be operated 
with ease as though the entire core and drum memory of the computer were indivi¬ 
dually addressable. The improved Group Control Program includes a modified group 
structure which is different from the old group structure which has been in use for 
some time. A Group consists of a number of Sections, each Section consisting of a 
collection of coding and data which performs some specific function. The new 
Group Control Program maintains a record of which Sections are in core memory and 
when a Group 1s called for, any Sections which are not already in core memory are 
brought in automatically for execution. Sections and Groups may be of arbitrary 
lengths and occupy arbitrary locations in the improved Group Control Program. 

The new program will result in greatly increased efficiency of computer application 
and increased flexibility as well. 

The Slurp program proper has been modified to Include several new facl 
lines; the Slurp program operates like a specially designed simulated computer 
with very unusual control features. Among these ure several modes of operations 
so that programs which differ only slightly may be Intermingled to occupy the same- 
storage space and yet perform different functions depending upon the setting of 
the control element of the simulated computer. The Slurp computer also operates 
Internally in terms of floating addresses of a special type, so that extreme flexl 
bllity in program modifications are possible. A basic feature of this Slurp program 
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is that its operation code includes sequences of arbitrary Flexo characters. 

When such sequences are encountered they are automatically uisplayed on any type 
of output equipment which has been previously selected and set up by other Slurp 
instructions. The control element of the simulated computer is arranged so that 
any word or sequence of words may be stored only once and yet arbitrary combinations 
of words and sequences can be assembled for display by making use of the Jump- 
type instructions which are built into the computer. In this way an extremely 
wide variety of labelling information can be produced with a minimum amount of 
computer storage. During the past quarter this feature has been used to provide 
labelling of results from special Slurp post mortem routines and scope plot 
routines entirely under the control of the Slurp computer. 

As described in the last quarterly report, an important feature of the 
Slurp system is a Logging Program which creates a logical record of all manual 
intervention actions and displays during a computer run on magnetic tape. During 
the past quarter an Editor Generator Program, originally coded by F. C. Helwlg 
of the S and EC Group, has been incorporated with Slurp. The Editor Generator 
Program represents another level of simulated computer, and its function is to write 
or generate editing programs which make use of the facilities of Slurp Itself to 
provide easily read edited versions of the information placed on magnetic tape by 
the Logging Program. T*o Editor Generator Program has an operation code which 
consists of individual F^exo characters so that It is attached to Slurp as though 
it were another output device. The editing programs themselves after they have 
been generated in tun. make use of Slurp to operate actual output equipment as 
above . 

D. T. Ross 

Servomechanisms Laboratory 

141 S & EC SUBROUTINE STUDY 

During the past quarter, a double register fixed point sign agreement 
routine has been tested and added to the subroutine library. This routine will 
test for the sign agreement of the major and minor parts of a number in 96 psec. 

If the signs disagree the sign of the minor will be made to agree with that of 
the major. If the major is equal to zero Its sign will be set equal to the 
sign of the minor. The routine takes a maximum of 512 pseconds. 

M. E. Callaghan 

Digital Computer Laboratory 

162 N. DETERMINATION OF PHASE SHIFTS FROM EXPERIMENTAL CROSS-SECTIONS 

In the last quarter we have been attempting to fit the phase shift com¬ 
binations in an area of our function containing an anomaly. We have found a few 
significant facts about this region leading to a decision to change our approach 
between the energies of 2.6 MEV and 2.72 MEV. Future plans Include using our 
present approach to determine fits beyond the energy of 2.72 MEV. 

Programmers working on this problem are B. Campbell and E. Mack. 

E. Mack 

Nuclear Science Laboratory 


193 L. EIGENVALUE PROBLEM FOR PROPAGATION OF E M. WAVES 

Calculations of eigenvalues, eigenfunctions and mode sums have been con¬ 
tinued lor the bilinear model at 410 Me and 3000 Me. Convergent sums have been 
obtained at four additional distances at 410 Me, and at the two largest distances 

used for 3000 Me. 

More accurate programs for obtaining eigenfunctions and eigenvalues were 
written for the inverse-square model. These calculations ure being made at 50 Me. 

R. M. Ring 
Lincoln Laboratory 

199 N. LAMINAR BOUNDARY LAYER OF A STF.AUY COMPRESSIBLF FLOW IN THE ENTRANCE 
REGION OF A TUBE 

Solutions were obtained for the case where temperature dependence of the 
viscosity and thermal conductivity of the compressible fluid is taken into con¬ 
sideration . 

More solutions are to be obtained for different entrance Mach numbers 
and thermal conditions at the tube wall. 

T. Y. Toong 
Mechanical Engineering 

203 C. RESPONSE OP A MULTI-STORY FRAME BUILDING UNDER DYNAMIC LOADING 

Time during the past quarter has been used in examining the effect of 
variation of certain parameters. A generalizing program to handle other struc¬ 
tures has also been worked on 

R. G Gray 
Civil Engineering 

219 TRANSPORTATION PROBLEM 

See report for problem 326. 

J. Dennis 

Electrical Engineering 

225 B.N. NEUTRON-DEUTEKON SCATTERING 

During the past quarter, calculations for a non-zero value of the polarl 
zatlon parameter, p e .02, with the expanded four parameter trial function, were 
made. As in the case of no polarization, these results differ little from the two 
parameter case. The next step will be to try a higher value of y. 

H. Paul 

Nuclear Science Laboratory 
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236 C. TRANSIENT RESPONSE Of AIRCRAFT STRUCTURES TO AERODYNAMIC HEATING 

For the past three months Whirlwind routines for solving secular equa¬ 
tions have been used to aid the solution of certain heat flow problems. 

L. A. Schmlt 

Aeroelastic and Structures 
Research Laboratory 

245 N. THEORY OF NEUTRON REACTIONS 

During this quarter, more tests were made in critical areas to deter 
mine better values of c and y . When these are determined a complete case will 
be done for the cross-sections and angular distributions. 

E. Campbell 
Nuclear Science 

253 N. APW AS APPLIED TU FACE- AND BODY-CENTERED IRON 

An investigation of the accuracy of the Noumerov method of integrating 
the radial wave equation has been completed. This check was carried out by in¬ 
vestigating the solutions of the Coulomb wave equation, for which the analytical 
solutions are known. 

Currently an empty lattice test for the face-centered structure is 
being carried out. Programmers working on this problem are M. M. Saffren and J. 
H. Wood. 

J. H. Wood 

Solid State and Molecular 
Theory Group 

256 C. WW1-1103 TRANSLATION PROGRAM 

IXirlng the past quarter the input language vocabulary for the transla¬ 
tion program has been extended by several modifications to the program. The form 
of octal and decimal numbers has been generalised to permit factors of powers of 
two and ten. The work FEED permits programmers to obtain enough blank tape 
where desired in the translated tape to meet the dimensional problems of the 1103 
input paper tape reader. 

Programmers working on this problem are F. C. Helwig and J. M. Franko- 

vich. 

J. M. Frankovich 
Lincoln Laboratory 

257 C. HORIZONTAL STABILIZER ANALYSIS 

Some production runs were made during the last quarter; others will be 
made during the next quarter. The program is also being used to investigate 
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several parametric variations. 

Programmers working on this problem are N. P. Hobbs and E. Crisclone. 

E. Crisclone 

Aeroelastic and Structures 
Research Laboratory 

261 C. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURES 

During the past quarter, the original problem was in the process of 
being finished. 

M. J. Buerger 
Geology and Geophysics 

273 N. ANALYSIS OF AIR SHOWER DATA 

During the past quarter about 1500 air showers have been analyzed by 
the WW1 computer. Very little programming time has been spent on the problem 
since the program works satisfactorily and has required only occasional minor 
modification. We plan to continue with the reduction of air shower data using 
the WW1 computer and also to develop a similar program for the new IBM 704. 

Programmers working on this problem are G. Clark and F. Sc herb. 

F. Scherb 

Nuclear Science Laboratory 

27 1 N. MULTIPLE SCATTERING OF WAVES FKOM A SPATIAL AllHAY OF SPHERICAL SCAiTERERS 

A satisfactory fit, F(k Z ,K) lias been obtained for the sum 2 
for all seven directions of K. F(k 2 ,K) takes the fora: 

(a D + b c K + c 0 K 2 + doK 3 + doK 4 ) + (ay + byK + cyK 2 + dyK 3 + eyK 4 ) k 2 
+ (a 2 + b 2 K + c 2 K 2 +■ d 2 X 3 + e^Mk 2 ) 2 + <a 3 + b 3 K + c 3 K 2 + d 3 K 3 + e 3 K 4 )(k 2 ) 3 
+ (a., + b 4 K + c 4 K 2 + d 4 K 3 + e 4 K 4 )(k 2 ) 4 . 

Work is now in progress to obtain fits for the remaining six sums. 

M. Karakashlan 
Physics Department 

278 N. ENERGY LEVELS OF DIATOMIC HYDRIDES Lii! 

In order to obtain a better approximation to tho true molecular wave 
function for the lithium hydride molecule, 2, v orbitals arc being included for 
both the lithium and hydrogen atoms. While radial correlation has been rensonably 
well considered in the treatment so far, no angular correlation has been included 
since Yf orbitals were omitted. 

The necessary integrals are being evaluated at the present Hint by 
means of the programs, written by Dr. Corbat^, which have previously been used. 
Several lnternuclear distances are being considered for which we have results of 
earlier calculations. The self-consistent field calculation will b carried 
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out using Or. Nesbet’i programs. If there Is time, a configuration Interaction 
treatment will be Included. 

A. M. Karo 

Solid State and Molecular 
Theory Group 

290 N. POLARIZABILITY EFFECTS IN ATOMS AND MOLECULES 

A large part of this quarter was used in analysing and organizing the 
results of previous computer works and in writing up the results as part of a 
thesis submitted to the M.S.T. Physics Ltepartoent. Most of the physical results 
and interpretation have been or will be reported in the Quarterly Progress Re¬ 
ports of the Solid State and Molecular Theory Group. 

We have recently started to use the computer again and an intensive pro¬ 
gram is planned for the remainder of S&J2C Whirlwind operations. There are several 
general points which can bo made in regard to the use of the computer on this pro¬ 
blem and it seems appropriate to give them here. Over the course of several years 
a large number of routines and sub-routines have been built up that are especially 
applicable to the solution of problems in the electronic structure of atoms, mole¬ 
cules, and solids. Most of these have been written by Corbato, Meckler and Nesbet 
(with additions and modlfications by others) and have been adequately described 
in Whirlwind I Quarterly Progress Reports. The capability of these programs have 
now been thoroughly explored and their use on a more or less production basis has 
resulted in significant advances in solid state and molecular theory. In parti¬ 
cular, in our works we have been concerned with ab initio electronic wave func¬ 
tions for atoms and simple molecules and these programs have enabled us to carry 
various aspects of the problem considerably further than has been possible before. 

In looking forward to the l.B.M. 704, initially we will be at a consi¬ 
derable disadvantage since competent programmers are not now available in suf¬ 
ficient numbers to transform these routines in a short time. It may be expec¬ 
ted that there will be a gap of one to two years in re-establishing certain of 
these programs. Because of this and because we have these routines in a form well 
tailored to our specific type of molecule problem, production runs In the next few 
months will be extremely valuable. It is reasonable to hope that this will allow 
us to obtain adequate theoretical predictions for experimentally inaccessible 
quantities of direct aid to spectroscopy and upper atmosphere physics. 

L. C Allen 

Solid State and Molecular 
Theory Group 

300 L. TROPOSPHERIC PROPAGATION 

A subroutine for generating incomplete elliptic integrals of the first 
and second kinds has been programmed and tested. There has been no other activity 
this quarter. 

W Mason 

Lincoln Laboratory 


JU9 B.N PUKE \NU IMPURE POTASSIUM CHLORIDE CRYSTAL 

Certain additional details of the structure of the valence band of KC1 
were obtained by solving seculur equations at new points in reciprocal space. 

Some calculations of the charge density in KCI have also been made. In the near 
future it ib planned to attempt a density of states calculation for the valence 
band of KCI and to obtain the charge density to a somewhat higher accuracy. 

L. P. Howland 

Solid State and Molecular 
Theory Group 

310 C. TRAJECTORY CALCULATIONS FOR A ROCKET DURING POWERED FLIGHT 

By eliminating intermediate printout, total computation time per run 
has been considerably reduced. Parametric variations are planned for the im¬ 
mediate future. 

L. J. Berman 

Aerophysics Research Group 

312 L. ERROR ANALYSIS 

A modification of 100-word length has been written and tested. Pro¬ 
duction runs are being made. There may be other modifications of the problem in 
the attempt to use the main program for satellite computations. Programmers 
working on this problem are E. Hutcheson and L. Peterson. 

L Peterson 

Lincoln Laboratory 

317 C. EXTRACTION OF STABILITY DERIVATIVES FROM FLIGHT TEST UATA 

In the past quarter, this project has been concerned primarily with the 
problem of extracting lateral stability derivatives. A number of simulated re¬ 
sponses to aileron and rudder inputs have been obtained. Generally these lesponses 
lie in the range within which aerodynamic non-linearities are unimportant, but 
for many of them, nonlinear inertial terms ( cross-coupling terms) are significant 
From these responses extractions of the derivatives have been made by both Shln- 
brot's method and a new statistical method developed by this project. The effects 
of varying the tyoe of response, the assumed precision of measurement, the number 
of frequencies used and other factors influencing the accuracy of results have 
been investigated. It is proposed to continue the present investigations in the 
next quarter and to obtain responses from which other quantities in the equations 
of motion of an aircraft which have not thus far been considered can be extracted. 
Among these additional quantities are l z , I xz . C Ba •, and parameters determining 
nonlinear lateral aerodynamic terms. 

Programmers working on this problem are M. E. Moult, L. L. Mazzola, M. 

N. Springer, L. E. Wilkie and J F. White. 

M. N. Springer 

Aerophysics Research Group 
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320 C. PRODUCTION FUR TilANSPOHTATIoN PROBLEM 

During the past quarter, additional data was obtained on Lite computation 
time used in solving li unsportation problems, when the best method of selecting the 
next step is used (s«*‘ -Summary Report No. 48). This data is presented in Fig. 1. 

No further computer use is contemplated for this problem. 



1 2 5 10 20 50 


Computation Time in Minutes 


Fig. 1 Computation Time Duta - Whirlwind Transportation Routine. 

m = No. of Supplies 
n - No. of Demands 


J. Dennis 

Electrical Engineering 

327 L. PREDICTION ANALYSIS 

The program has been completed and production runs are being made. Modi 
ficatlons may be made from lime to time. Programmers working on this problem are 
L. Peterson, E. Hutcheson and I*. Willmann. 

E. Hutcheson 
Lincoln Laboratory 

337 N. NONLINEAR SECOND ORDER DIFFERENTIAL EQUATIONS IN THE THEORY OF ELASTIC 
SHELLS 

The |>rottram evaluHting Physical quantities of Interest In the present 


WMIKLWIND CODING AND Al I'l.K AlloNS 

problem, which wut> mentioned in the last progress report, has been successful I > 
used to obtain stresses und displacements for a wide range of parameters of the 
differential equations in question These results are included in a report written 
up i' r publication. Some experiments have also been made during the last quarter 
with .i moditied program in order to obtain "asymptotic solutions of the differen 
r»a, equations subject t»» suitable boundary conditions. However, no reliable re 
suits have been obtained as yet. We feel that further modifications of the pro 
sent program will be necessary to get the desired results, a point unuer invest! 
gat Ion at the present time. 

H. Weinitschke 
Mathematics Department 

341 G. STATISTICAL AND DYNAMIC METHODS IN FORECASTING 

Studies of sea'level pressure by empirical orthogonal functions con¬ 
tinued. Certain functions of pressure and pressure tendency, selected by a sta 
tistical screening procedure, were used as predictors in a linear forecast method. 
Programs were developed for Whirlwind to test this new method on a set of inde¬ 
pendent data. 

1 method for predicting the surface 5-day mean temperature anomaly based 
on the coefficients of the empirical orthogonal functions was tested on an inde¬ 
pendent sample. Verifications of the resulting forecasts are now in progress. 

Work continues on a program for a simplified numerical model of the 
general circulation of the earth's atmosphere developed in terms of spherical 
harmonics. 

Another phase of the Statistical Forecasting Project has been devoted to 
the development of a stat1st leal-dynamic method of numerical weather prediction 
Work on Whirlwind has been concerned mainly with the determination of certain con 
stunts appearing in the prognostic equations. 

Ail of the work done on the project has been supervised by Professor 
Edward N. Lorenz, Department of Metoorology. Programmers working on this problem 
are E, A. Kelley, U. Shorr and K. Bryan. 

E. A. Kelley 

Metenrol ogy Department 

343 C. WEATHER PREDICTION 

During the past quarter, the computations on this problem were Completed. 
I he results of this work way be found in the final report Resear ch o n At nos)uci1«- 
Pressure Changes . 

E. A Kelley 
Meteorology Depulimont 

:;.'i0 D. COMPUTATION OF VARIANCES AND COVARIANCES 

Empirical functions of temperature hnve been computed and their !»•:• and 
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contemporary correlations with empirical functions of pressure found. No further 
machine computations are planned for the immediate future. 

Programmers working on this problem were B. Shorr and D Gilman. 

D. L. Gilman 
Meteorology Department 

354 D. RESPONSE OP A SINGLE STORY CONCRETE BUILDING 

Airing the past quarter the problem was completed. 

S. Namyet 
Civil Engineering 

361 B,N. GROWTH OF FATIGUE CRACKS 

Airing the past quarter, the integral equation expressing crack growth 
was debugged and several runs were made. A program to determine the stress die 
trlbutlon In a round bar with a circumferential crack in It was set up and is 
now being debugged. 

J. B. Walsh 
Mechanical Engineering 

364 C. BLAST RESPONSE OF ROTOR BLADES 

Airing the past quarter routines were debugged. In the future it Is 
hoped that the equations of motion will be integrated. 

A. Holland 

Aeroelastics and Structures 
Research Laboratory 

371 L. ATMOSPHERIC PROPAGATION OF RADIO WAVES 


WHIRLWIND CODING AND APPLICATIONS 

377 L. COVERAGE ANALYSIS 

A scries of production runs la being made. Additional modlf1cat Ions 
may be made In the future. Programmers working on this problem are L. Peterson 
and K. Hu t cheson. 

E. Hutcheson 
Lincoln Laboratory 

380 B. SWITCHING CIRCUITS 

A program for simplification of switching functions of 14 or fewer vari¬ 
ables by means of the McCluskey-Quine reduction procedure has been completed and 
checked out and 89 problems have beon solved. The switching function to be sim¬ 
plified is speciflod In standard sum form The program compares the torms In 
the standard sum to form the prime impllcanta and then selects certain of these 
prime lmplicants to form the minimum sum 

A program has been written for synthesis of nor-serles-parallel contact 
networks by using Boolean matrix methods. The standard sura which represents the 
transmission of the network to be synthesized Is represented by an octal fraction 
By means of successive node insertions, a Boolean matrix is obtained which re¬ 
presents a network having the desired transmission. After elimination of re¬ 
dundant contacts, the matrix representation of the contact network is printed 
out. In approximately half of the 67 synthesis problems which have been solved, 
minimum networks were obtained. In most of the remaining cases, the networks con¬ 
tained only one or two extra contacts. The program le capable of synthesizing a 
network to represent any four-variable switching function, and the solution time 
ranges from a fraction of a second to several minutes. 

The computer work on this problem is practically completed. After 
solution of a few more problems with each program, the results will bo written 
up In a Master's thesis entitled, Digital-Computer Solution of Switching- 
Circuit Problems. 


A rapid method of reasonable accuracy for integrating the range-height 
equation (described In earlier guarterly Progress Reports) has been programmed. 
T Is program will be used with large sets of experimental data and statistical 
studies Will be carried out on the results. 


W. Mason 

Lincoln Laboratory 

372 B. DESIGN OP SPHERICAL SHELL SEGMENTS 


C. Roth, Jr. 

Electrical Engineering 

38G C. FREE CONVECTION 

Additional similarity solutions have been completed corresponding to 
several choices for the constant* defining the surface temperature variation 
(~BX^, X being the length parameter in the flow direction). Final computations 
are expected to be complete within the next quarter. 


During the past quarter final results were obtained and the full purpos 
of the program was accomplished. No future work on this program is necessary. 


£. Traum 

Civil and Sanitary Engineering 
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J. R. Baron 

Naval Supersonic Laboratory 


387 C. DETERMINATION OF VELOCITY POTENTIAL 

Tho complete program for the delta wing has boen debugged and 1* 
available and has been used for production runs to yield useful Information. In 
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addition, the program has boon expanded in scopo so that the machine takes the 
basic aerodynamic influence coefficients and combines them with the proper weighting 
(actors to form the various component matrices used in the deterrumut ion of the ve 
ioclty potentials. Previously this was done on dusk calculators and was a tedious, 
lme-consuming, and error-prone process. The increase in lime per production run 
for inclusion of this stop is a small fraction of the total machine time required. 

At present, three reduced frequencies at the lowest Much number and one reduced 
frequency at the intermediate Mach number have been investigated for the delta. 
Remaining are runs for two reduced frequencies at the intermediate M and 3 re 
duc<d frequencies at the high M. 

A first draft of the complete program for the straight wing has been 
completed and is in the process of being debugged. The swept wing program a 
waits the success!ul conclusion of the delta and straight wing programs. The 
swept und straight programs are essentially merely modifications of the succes¬ 
sful delta program. For each of the threw plunforms, three reduced frequencies 
must be investigated for each of three Much numbers. 

H. Staple!ord 
S. Gravity 

Aeroelastic and Structures 
Research Laboratory 

388 D. TEMPERATURE DISTRIBUTION IN AIRCRAFT GENERATORS 

A routine which produces a generator rating chart from design parameters 
has been written and is being tested. This program, when checked out, will be 
the final result of the present problem. 

R. M. Moroney 

Servomechanisms Laboratory 

389 D. SUPERSONIC FLOW OF AIR IN A TUBE 


Problem Description 


Physical Model . A supersonic stream of air enters a round tube with a uniform 
Mach number across the stream. Gas, which may be different from that in the 
entering stream, is forced through the wail which has uniform porosity. The 
flow in the entrance length of the tube is assumed to consist of two concentric 
regions: a boundary layer and an isentropic core. 

Kequlred Solution . Profiles of velocity, temperature and injected gas concentra¬ 
tion are required as u function of the distance from the entrance of the tube. 
Pressure distributions, friction factors and heat-transfer coefficients are also 
lequlred. The effect of entrance Mach number, foreign gas molecular weight, in¬ 
jection rate and heat transfer rate on the above quantities is required. 

Basic Differential Equations . The partial differential equations describing the 
problem are obtained from the Navier Stokes equations tn cylindrical coordinates 
by an order of magnitude analysis. These equations are as follows: 
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Overall Continuity 


_B («y>u) + B_ (r/> w) » 0 


Momentum or Eauatlon of Motion 


fu_» + w Url =1 1 &_ fru 3») - 6p 

' l hr '■* | Re r Sr |_ hr / 


Single Species Continuity 


8e 0 a f u Be + w be 1 = 1 S_ f r>> r&c + 

2 \ hr bz J r hr l [hr 


ac(l-c) b(lnT) 

6 r J 


Energy 


KeonCp f u 6T + » St] = Re 0 M§ (y 0 -l)w + 2_ 1 5_(r LhT ) 

l 5r bz J fjz Pr 0 r &r } hr J 

+ ?_ 2 fi* + OWW-O) 1 §li &_ f T [Cpi-Cp 2 + 7^,-1 0*21 |] 

Sc Q r L&* T 6t| l 7o c(M 21 -1> +ljj 

+ To-1 *2lT He n j»Cu he + w be") + Mg (y 0 -l) p (b w\ 2 
y 0 e(M 21 -ll +1 L hr hz | \br) 

5. Equation of State and Properties 

P = ^oRT 
h = CpT 

In these equations all lengths are referred to the tube radius, a,all 
velocities to the initial free stream velocity, w Q , all properties to free stream 
entrance properties, uor ^i2 0 • and the pressure to The subscript 0 re¬ 
fers to entrance free stream conditions. 

Tran sformat ions . The overall continuity equation is satisfied by the introduction 
of a stream function: 

/our = - Ft 


A modified length Reynolds Number rand a wall distance parameter 
are defined by: * 


7 X*) . 

7 ' _; i ['' ( : ) 1 


Re 0 = 2wa 
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The boundary layer quantities * , T, ft, C are expanded in aeries in J 
with coefficient which are functions of 

* - - f i t j ) 

t « t q ii 1_ «i (•y) 


c ’ fi s 1 ’ 1 gi< ? ) 


The properties in the lsentroplc core are described by series injr . 

W c - W Q (1 + HZ r l K i> 

ini 

P = Po <1 + ^ J 1 
i=l 

T ' T (1 + 2C / t A ) 

1=1 

p » />(1+ SI Y 1 

' 1=1 J 

The coefficients f or p , T c , a,/> c ““y b ® related to the coefficient by 
use of the lsentroplc relations for the core. 

When the series for T, <f , ft, C are substituted into the partial dif¬ 
ferential equations, the coefficients of the powers of J in each equation may 
be equated to aero, since these equations must hold anywhere in the region of 
validity of the analysis. Three infinite sets of ordinary differential equations 
in fy. gj, ft if are obtained. Each set of three determines the coefficients 
of one power of J in the series. The equations for the first terms are: 

*1 ' - - t pl + 2 *i’ - f e i' + °i ti 

l M e l J V hi j 9 1 

«1 = - (Pr 0 l i + A{) OJ_ - Pr„ m (r 0 ' IHfjOj + fi"®!) 2 

\. -1 K< 

where: ^ = jOljT + 2 T9 + 3 T0) , = j 0"( 1 T + 2 T9 + 3 T0) 

H- = jO^T + 3T»)9' 3 = i<^( 2 T + T 3 0)0' 

and , ^er, 1 T, T, § 3 T are all constants. 

Only the three first order equations are nonlinear. The higher order 
equations are linear and contain functions calculated from the lower order equa¬ 
tions, but their complexity increases rapidly with the order of the equation. 

The boundary conditions for the first order equations are: 

•y « 0; f t « 0, fj = 0, • 0, = 0 

vis oo ; limf* = 2, lin 8 * 1, ling = 0. 

* *« 1 7n 1 1 
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Nu merical Solution . The numerical part of the solution to be computed on Whirl 
wind will use a subroutine to calculate the highest derivatives in each of the 
equations for arbitrary values of the variables in conjunction with the Kutta 
Gill integration method. Unknown values of 6}, gj, fi f'/ must be guessed at a 
value of y equal to zero so that f^, , fi gj will approach the proper limits 

at large values of^ . Two initial calculations are made, then linear interpola¬ 
tion is used to establish the conditions for the third and successive iterations 
based on the two most recent values. When the boundary conditions at the Be 
cond point are satisfied within the prescribed tolerance, the quantities which 
depend on the profiles such as heat-transfer coefficients and friction factors 
are calculated. 

The linear equations for the higher order terms are solved by obtain¬ 
ing particular solutions which satisfy the conditions at the initial point and 
adding to these solutions which contain parameters which may be adjusted to 
satisfy the second point conditions. 


A routine has been programmed and checked to calculate the first 
terms in the series solution. This method is dependent for convergence on the 
trial values of initial conditions with which it is started. Therefore, ex 
ploratory calculations or a modification of the method will be necessary. 

The routine to calculate the second terms in the series Is now being 
programmed and a hand calculation being made to check it. Calculation of a 
third set of terms may bo undertaken if required for convergence of the series 

J Radblll 

Mechanical Engineering 

390 Q. MITCHELL S WAVE-MAKING INTEGRAL 


Problem Description 

Mitchells integral expression in the form 


i =/« 


cos(xt) t dt 


expresses the mutual wave resistance between any two elements on the hull of 
a ship (progressing at s uniform velocity) when multiplied by tne mean areas 
and the mean slopes of the two elements concerned. In the above expression: 

t Is the variable of integration 

y is a function of the vertical spacing of the areas and of speed 
x is a function of the horizontal spacing of the areas and of speed. 

At the lower limit this integral is infinite as are also all previous trans¬ 
formations of this Integral. No analytical expression for this integral has 
yet been found However, a transformation has been made of this Integral to 
the form: 
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wo 

1 = 1 /e~ t/y fj. (2t" + x) + J. <2t '/*■ - x)7 dt 

l 1 i 

in which is the Bessel function of the first kind of order one. This trans¬ 
formation is finite at all points (y is never zero for a finite area), thus enabling 
numerical integration to bo performed without the use of questionable formulae 
near the discontinuities. The convergence of the integrand is reasonably rapid 
in all cases for such an evaluation to be practical. With a constant x and cons tan 
y, t can be varied in small increments and use of Simpson's rule then gives the 
value of the Integral. This is the part of the problem that Whirlwind can perform. 
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Reports 33-43) a library of subroutines was developed to assist in data preparu 
tion for the MIT NCMM Under problem 250 (Whirlwind Progress Reports 41, 42) the 
feasibility of applying modern automatic programming techniques to machine tool 
programming was successfully demonstrated for plane parts bounded by straight line 
and circular segments. Problem 394 1 b concerned with the extension und cluboru 
tion of these automatic programming techniques to arbitrary surfaces which can be 
machined by 3 and 5 hxIb numerically controlled milling machines. The programs 
are currently being developed by Douglas T. Ross, Harry E. Pople, Jerome K. Wonker, 
George Zames. Samuel M. Matsu, Servomechanisms Laboratory Staff members. This 
work is sponsored by the Air Materiel Command Through DSR Project 6873. 


Progress 


The final form used for numerical integration with Simpson s three 
point rule is 


- V 


> <l2+1 > 2 c OS x(. 2 + 1 > (t 2 + l> = 


dt 


J(t 2 + 1) + l 


The problem may be viewed as the design of an APT System (Automatically 
Programmed Tool system) combining the human operator, the machine tool itself, and 
an APT computer, which is a combination of a general purpose computer and existing 
machine tool directors. The objective is a system of Whirlwind computer programs 
which will simulate an APT computer with a language which is convenient for hu¬ 
man expression of machined part designs A design expressed in this language 
will be automatically translated by the APT computer into machine tool instructiorts 


rather than the initially proposed transformation containing Bessel functions. 

Values of the above integral expression have been obtained on WWI 
These values apply to ship forms that are less than 400 feet and advancing at 
normal speeds. The remaining work is centered around a refined grid in the 
neighbourhood of the small region whore the integral takes on large values and 
obtaining values of the integral for ships longer than 400 feet. 

This data is being incorporated in Table II of my Sc.D. thesis en¬ 
titled The Slope-Area Concept fur determining the Wave Resistance of a Ship . 

W. Tudor 

Civil Engineering 

391 L. MAGNETIC RELAXATION IN THIN FILMS 

The programming has been perfected so that future calculation only in¬ 
volves parameter change (also, three runs have been made bearing on current ex¬ 
periments). The amount and exact form of future calculation will depend on the 
course of the concurrent experimental study of film relaxation. 

Programmers working on this problem are J. Frankovich and D.0. Smith. 

0. O. Smith 

Lincoln Laboratory 

394 C AUTOMATIC PROGRAMMING FOR NUMERICALLY CONTROLLED MACHINE TOOLS 

T\\e completion of the MIT Numerically Controlled Milling Machine in 
1952 introduced a new ora in the automatic translation of design requirements 
Into finished parts. The primary problem associated with this type of machine 
tool, the preparation of detailed instructions for making a pari, has been stu¬ 
died using the Whirlwind Computer. Under problem 132 (Whirlwind Progress 


During the past quarter coding and debugging has continued on Whirlwind 
programs which allow programming of a part by means of arbitrary curves on the 
surface of a part, and also in terms of entire regions of the part. A number of 
additional functions have been put into the system programmed by curves, such ns 
automatic indexing of curves to allow a specialized type of region programming. 
The oscilloscope output routine which draws near perspective views of the part 
being made has been incorporated into these routines and a number of test parts 
have been drawn automatically by the programs. Progress on Problem 394 was de¬ 
monstrated in March as part of a special course on Programming f r Numerical Con 
trol presented by the Servomechanisms Laboratory staff under the auspices of the 
Summer Session Office. Work is now underway on more general programs. The de¬ 
sign of special purpose, English-like languages for programming A1‘T Systems will 
be begun in the next quarter. 

Associated studies are being made for the determination of appropriate 
languages for the feedback of information from the APT computer to the human op 
erator. A progrum has been written by which the Whirlwind Computer can draw on 
the output oscilloscope an axonometric projection (near perspective) of the part 
being made. The program allows arbitrary orientation of axes and arbitrary 
scaling so that any view of any portion of the part can be shown The program 
has several modes of operation, but at present the picture is displayed us the 
path of cut vectors swept out by the motion of the cutting tool over the surface 
of the part. English-like computer to human languages will also be devised. 

D. T Ross 

Servomechanisms Laboratory 

397 N RESPONSE FUNCTION OF AIR SHOWER DETECTORS 

The response of the detector array used on the roof of Building 6 has 
been completed. An arrny now in use in India has been analyzed in part. An ar 
ray for use in the study of altitude variations of air showers will be analyzed. 

0. M Clark 

Physics Department 
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399 L.B. DOMAIN WALL MOTION 


Proble m Des cription 

In studying the motion of 180° magnetic domain walls In magnetization re¬ 
versal processes In ferrite toroids an equation of the form: 

x + 2d x + f(x) * h 

occurs where x is the normalised distance the wall moves, 2d Is a damping parameter 
f(x) is a nonlinear restoring function with maximum and minimum amplitudes of *| 
and is periodic in x with period 2, h is a normalized driving force. The problem 
Is to obtain solutions of this equation for various values of d and h. In par¬ 
ticular a plot of x or wall velocity vs. time is desired. It Is expected that 
such plots will have the same general character as the actual voltage response 
observed in magnetic cores during switch processes under the Influence of vary¬ 
ing driving currents. 


Cp specific heat. 

T temperature of adjoining elements 

b,c,etc * 


St ress Equation s: 

m m 

°*i B i C i* T l + Vie S E l L l* T i X ie A l 
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Progress 

Ujrlng the past quarter, the initial computations were completed. It 
is now planned to reprogram the problem along somewhat different lines. Pro¬ 
grammers working on this problem were J. V. Harrington and G. Dineen. 

J. V. Harrington 
Lincoln Laboratory 

400 C. TEMPERATURE AND STRESS RESPONSE 


The number of repetitions of the temperature equations depends on the 
time increment chosen and the time values of interest. The temperatures are 
stored at the times of interest for use in the stress equations. The stress e- 
quations are repeated only for the time points of interest. 

Progress 

During the past quarter, the program has been modified to include the 
response to a thermal pulse. Programmers for the problem are J. C. Lorla and J. 
F Duvivier. 


Problem Description 

This problem is concerned with the temperature and stress responses of 
various helicopter blades exposed to thermal radiation on one surface. The ther¬ 
mal inputs include a step input and time-varying pulses. The mathematical ana- 
lysis consists of a second order partial differential heat conduction equation 
which describes the absorption and diffusion of heat energy and the Timoshenko- 
Goodler stress equations for predicting the thermal stresses Involved. These 
equations, reduced to finite difference form for solution by the Whirlwind com¬ 
puter are as follows: 


Temperature tjq untlons : 


r i.j * * Al< A l n w > + ElTj.j-i + OiTbj.x + HjTc.j-i + FiTd,.! + IyTe.j-1 

lcr- c p 

where K l + Gl + H 1 4-F i -f I A =1 

and T i , J = temperature in element 1 at time J. 


A 'itj< A texp> = heat absorbed by element i in the time interval t. . 

where A iex p is the exposed surface area of element 
Al cr* cross-sectional area of element 1. 
y — weight density of element. 


t° tj 


J. C. Lorla 

Aeroelastics and Structures 
Research Laboratory 


401 N. NON-STATIONARY QUEUEING PROBLEM 


Problem Description 

1. Let Pi feio **H • • -J • * = I. - . .1400/h be a set of discrete dis¬ 
tribution functions denned recursively by 

P, » *3 l> .8, for J mod 1440/h 

in J=0 1-1, n+l-J J.l 

where S = Q.jh)^ e~M h , where the ky and h are specif ted. the \y being periodic 
J.l j! 


The problem is to calculate the stationary vectors P*. 

. th 


with period i * ***0 
h 

2. The application is : 1) k± is the arrival rate in the l fc “ interval 
of a calendar day of airplanes at an airport; 2) h is the servicing time of an 
airplane and is either 1 or 2 minutes as specified; 3) l'i n is the probability 
that at the end of the l ttv interval there are n-1 airplanes waiting to land and 
one landing, unless n»0 in which case none ore waiting or landing. 

3. The Sj are computed in (24,G), the P n in (15,0). 
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Progress 

During the past quarter the problem was successfully completed. Eight 
production runs were made. 

H P. Galllher 
Operations Research Project 

402 N. MONTE CARLO INVENTORY CONTROL STUDY 

Problem Description 

The problem includes various studies of the effect of variability of 
demand for a stored product and variability of replenishment times, on the in¬ 
ventory and ordering frequency. The studies Include both single and multiple 
warehouse problems, the latter being an actual case study in the Ordnance Corps, 
U.S. Army, the former theoretical studies. 

Random numbers are generated by a multiplicative congruential tech¬ 
nique. All probability distributions are sums or probability combinations of 
independent exponential variates. All computation is done in (15,0) except for 
natural logarithm computations 

Progres s 

IXiring the past quarter two programs were written. It is planned to 
start trouble shooting in the immediate future. Programmers working on the pro¬ 
blem are H. P. Galllher and M. Simond. 

H. P. Galllher 

Operations Research Project 

104 D. CORE OPTIMIZATION 

This pioblem is an extension of Problem 270 with the following changes: 

1. Core materials to be low enrichment uranium and water moderator. 

2. Cladding materials stainless steel, zirconium and aluminum. 

3. Water reflector. 

4. Cylindrical geometry only. 

5. Experimental data from Brookhaven and Bettis critical experiments and 
EBWR. 

(See Quarterly Progress Report No. 16 tor a description of Problem 270.) 

Progre ss 

Time during the past quarter was used in the checking of results against 
experimental data and adjusting questionable parameters. Programmers working on 
this problem are R. Gardner and B. A. bee. 

B. A. Lee 

Nuclear Engineering 


APPENDIX 

1. VISITORS 

Tours of Whirlwind 1 installation Include a showing of the film, Making 
Electrons Count , a computer demonstration, and an informal discussion of the ma - 


Jor 

206 

computer components 
people visited the 

During the past quarter, the following 6 
Laboratory. 

groups totalling 


January 9 

MIT Course. Machine-Aided Analysis 

48 


January 15 

MIT Course. Introduction to Digital 
Computers - Coding and Logic" 

11 


January 28 

Lexington High School 

14 


March 4 

Belmont Hill School 

13 


March 11 

Librarians from Harvard University 
and Radcllffe College 

57 


March 28 

MIT 10-day Course, Programming for 

Numerically Controlled Machine Tools 

63 


The procedure of holding Open House at the Laboratory on the first Tues¬ 
day of each month has continued. A total of 38 people attended the 3 Open House 
tours during the quarter, representing members and friends of the MIT community, 
brown University. Remington Rand, Westlnghouse, Philips-Endhoven in Holland, and 
Boston Latin School. 

During the past quarter, there were also 51 individuals who made brief 
tours of the computer installation at different times. Represented by these indi¬ 
viduals were the Mathematical Center in Amsterdam, Georgia Institute of Technology, 
United States Air Force, Royal Canadian Air Force, Lockheed, National Cash Register 
Company, Max Planck Institute in Germany. Norden-Ketay, Amherst High School In Buf¬ 
falo, Stanford Research Institute, Stone and Webster, Radio Corporation of America, 
Christian Science Monitor. RIAS, Houghton Mifflin Publishing Company, General Radio 
Company Yokogawa Electric Works and Japan Overseas Company. 


2. ACADEMIC 


Introduction 

During the spring term 1957 there were a number of graduate subjects in 
automatic computation, numerical analysis, and electronic data processing offered at 
MIT. In particular, the following subjects related directly to machine computation 
were offered during the spring term: 


28 


29 
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Subject 

Description 

Units 



Year 

Instructor 




6.535 

Digital Computer Coding and 


3 


6 


4 

Arden 

VEltSONNEl. ON THE 

PROJECT 



Logic 











6.538 

Electronic Computational 

3 

- 

1 

- 

5 

G 

Verzuh 

COMMITTE ON MACHINE METHODS OF COMPUTATION AND 

NUMERICAL ANALYSIS 


Laboratory 











6.54 

Pulsed-Data Systems 


3 

- 

6 


G 

Widrow 

Faculty Supervisors: 



6.568 

Switching Circuits 

3 


1 

- 

6 

G 

Caldwell 




15.542 

Management Information Systems 


3 

- 

6 


G 

Gregory 

Philip M. Morse, Chairman 

Physics 


M3 5 2 

Advanced Calculus for Engineers 


3 

- 

6 


G 

Hildebrand 

Samuel 11. Caldwell 

Electrical 

Engineering 

M3 9 

Methods of Applied Mathematics 


3 

- 

9 


G 

Hildebrand 

Stephen H. Crandall 

Mechanical 

Engineering 

M412 

Numerical Analysis 

3 

- 

2 

- 

7 

G 

Hildebrand 

David Durand 

Industrial 

Management 

N22 

Nuclear Reactor Theory II 


3 

* 

9 


G 

Clark 

Herman Feshbach 

Physics 



It is apparent Iron the above list that these subjects are predominantly Graduate A 
subjects - indicating that the subject is primarily for graduate students. However, 
a large number of undergraduate students in their Junior and senior year do take 
these subjects on an elective basis. 

Electronic Computational Laboratory - 6.538 


Jay W. Forrester 
Francis B. Hildebrand 
Chia-Chiao Lin 
J. Francis iteintjes 
Ex-Officio 


Industrial Management 

Mathematics 

Mathematics 

Electrical Engineering 


The subject matter offered in 6.538 during the spring term 1957 was 
altered to provide the student with more time on the 650 Magnetic Drum Calculator. 
Specifically, each student had £U leas t 8 hours of time on the 650 machine. In 
this manner a number of very proficient programmers were developed by the end of the 
term. A complete description of the course content including problems, term papers, 
and other assignments will be provided in the next quarterly report covering the 
period April - June. 

The material previously offered on the Card Programmed Calculator has been 
discontinued and instead this portion of the course has been devoted to a study of 
the physical and functional features of the Type 704 Electronic Data Processing 
Machine. Specifically, the students were taught a limited amount of programming on 
the 704 machine. Unfortunately, however, the machine was not in operation in time 
to allow them to use the machine to carry out actual machine runs prior to the com¬ 
pletion of the term. Again, a detailed description of the 704 material, physical 
plant, etc., will be given in the succeeding quarterly report. 

F. M. Yerzuh 
Assistant Director 
MIT Computation Center 


Electrical Engineering 
Computation Center 
Computation Center 

PROJECT WHIRLWIND 

Frank M Verzuh, Head 
Dean N. Arden 
Sheldon P. Best 
Marion Callaghan 
Frank C. Melwig 
Jack Koseman 
Arnold Siegel 
Murray Watkins 
Munroe R Weinstein 


Dean N. Arden 
Fernando J . Corbato* 
Frank M. Verzuh 
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